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Abstract The measurement of electron density is important for medical diagnosis and charged particle radiotherapy 
treatment planning. Traditionally, electron density is obtained by CT imaging using the relationship between 
CT-number and electron densities established beforehand. However, the measurement is not accurate due to the beam 
hardening effect. In this paper, we propose a simple and practical electron density acquisition method based on 
dual-energy CT technique. For each sample, the CT imaging is conducted using two selected X-ray energy from 
synchrotron radiation. A post-processing dual-energy reconstruction method is used. Linear attenuation coefficients of 
the scanned samples are obtained by FBP reconstruction. The effective atomic number and electron density are got by 
solving the dual-energy simultaneous equations. Different phantoms and breast tissues were scanned in this 
experimental study under 10 keV and 30 keV monochromatic X-rays. The distribution of effective atomic numbers 
and electron densities of the scanned phantoms were obtained by Dual-energy CT image reconstruction, which agrees 
well with the theoretical values. Compared with conventional methods, the measurement accuracy is greatly improved, 
and the measurement error is reduced to about 1%. This experimental study demonstrates that DECT imaging based 
on synchrotron radiation source is applicable to medical diagnosis for quantitative measurement with high accuracy. 
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1 Introduction 


Nowadays, computed tomography (CT) has become a 
mainstay for medical imaging and nondestructive 
testing in modern radiology. The internal structure of 
an object can be exactly reconstructed from CT 
scanning projections. The electron density is important 
in photon and charged particle radiotherapy"), 
Traditionally, this information is obtained by CT 
imaging using the relationship between CT-number (in 
Housfield units) and electron densities established 
beforehand”!, However, a conventional CT system, 
which uses X-ray tube, inherently has a broad energy 
spectrum. The resulting CT number has some 
uncertainties due to beam hardening effect. Thus, the 
accuracy of electron density measurement is degraded. 

The dual-energy CT (DECT) was proposed by 
Alvarez’! in 1976. Using X-ray under two different 
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tube energies, the information of electron density and 
effective atomic number of inspected objects can be 
derived. This technique has been widely used in 
medical and industrial applications, especially in 
safety check and material discrimination. All the 
proposed dual-energy reconstruction methods can be 
broadly classified into three groups: pre-processing, 
post-processing and iterative reconstruction methods!” 
In polychromatic X-ray DECT systems, the pre- 
processing methods are widely used to avoid beam- 
hardening effects. The post-processing reconstruction 
methods have their advantages in computational 
simplicity. The iterative reconstruction methods based 
on statistical model and non-statistical model yield a 


high signal noise ratio!. However, the great 
computation cost decreases their practicability. In 
China, the different groups focus on the DECT 


techniques. Zhang et al. [6] proposed a practical DECT 
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reconstruction method based on spectrum estimation, 
thus obtaining the distribution of effective atomic 
numbers and electron densities from polychromatic 
X-ray projections under two tube voltages. Most of the 
research groups focused on the clinical application and 
evaluation of DECT". 

Because X-rays from synchrotron radiation 
source are highly monochromized, the effects of 
beam-hardening in CT imaging can be ignored. Hence, 
synchrotron radiation CT (SRCT) can better detect 
low-contrast lesions and obtain better image 
1997, 
Dilmanian et al.®! conducted single and multi-energy 
CT at the National 
Brookhaven, USA. They compared the image noise 
between SRCT and conventional CT. In 2010, Chen et 


al."! measured the linear attenuation coefficients of 


quantification than conventional CT. In 


Synchrotron Light Source, 


surgical breast tissues by SRCT in the energy range of 
15.0 to 26.5 keV at the SYRMEP beamline of the 
ELETTRA SR facility in Trieste, Italy. Torikoshi et 
al." °"! and Tsunoo et al."'*! presented a method for 
DECT using monochromatic X-rays from Spring-8 in 
Japan, thus claiming to measure the electron density 
with accuracy of around 1%"). In China, the SSRF 
(Shanghai Synchrotron Radiation Facility) is a third- 
generation of synchrotron radiation light source. With 
the electron beam energy of 3.5 GeV, it is the fourth 
energy high in the world and offers X-ray with high 
intensity, high brilliant, high polarization and high 
stability!’ and it serves as an scientific platform for 
science research and technology development! 6", 

In this 


measurement method is proposed. A DECT imaging 


paper, a simple but practical 
experiment is designed using synchrotron radiation to 
explore its potential medical applications, and carried 


out at beamline BL13W 1 of SSRF. 


2 DECT reconstruction method 


Synchrotron radiation source provides highly 
monochromatic X-rays with continuously adjustable 
energies. Thus, the pre- and post-processing dual- 
energy reconstruction methods give similar results. For 
the purpose of simplicity in implementation, a post- 
processing dual-energy reconstruction method was 
used in our experiment. First, linear attenuation 


coefficients of the scanned samples are obtained using 


CT reconstruction using the two selected X-ray 
energies. Then the effective atomic number and 
electron density are solved from dual-energy 
simultaneous equations. 

In the photon energy range normally used in 
medical imaging, the linear attenuation coefficient u 
can be linear 


approximately expressed as a 


combination of photoelectric interaction term and 


Compton scattering term as follows" ®, 


U~O0.5K Naf E)Z" pe+0.5K:N4fk(EJpe (1) 


where f(E) and fxn(E) as functions of the X-ray 


energy are the energy-dependent photoelectric 
interaction and Compton scattering, respectively. The 
n (=~4) is a constant. The variable Z is the atomic 
number, and Na is the Avogadro’s number. K; and K3 
are constants related to the photoelectric interaction 
and the Compton scattering. pe is electron density 
(g/cm°) and can be expressed as 


Pe =2pZIA (2) 


where p is the mass density (g/cm?) and A is the 
atomic weight. For compounds, the effective atomic 
number (Zerr) is used instead of Z, and defined as 


I/n 
Log = Dez" 1X cZ, | 
l (3) 


where c; is the number of atoms for the i” element in 
the molecule of the compound and Z; is the atomic 
number of the i” element. Eq.(1) can be simplified to 


LEJ afl E)Z pe +02 fxn(E)pe (4) 


a, and az as constants are defined as follows. 


aı=0.5KıN4, a2=0.5K2N4 (5) 

Scanning with dual X-ray energies E; and Ep, 

the linear attenuation coefficients u(Ez) and u(Ep) can 

be obtained by the filtered back-projection (FBP) 

reconstruction method. We can derive dual-energy 
simultaneous equations from Eq.(4), as follows 


ME1)= af lEL)Z pe +02 fxn(E1)pe 
ME n)= afa Enz" Pe +02 fxn(Ex)pe (6) 


Solving the simultaneous Eq.(6) gives Eq.(7). 


Zett “laz fxn(Ex) MEn)—@2 fxn(En) MED) [ai fon(Ex) 
u(Er)} a Son(Ex) MEn)| (7) 
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P=[MEr) WE taal fon(Ex) Spn(En) Zee” +02[fxn(Ex) 
Sxn(En)]} (8) 


This method is available for low Z elements, 
but not high Z elements because they have absorption 
edges in the used energy region. Simulation study has 
been conducted by a numerical phantom consisting of 
different materials, indicating that the electron density 
can be obtained by high accuracy using the proposed 


dual-energy reconstruction method!'?!. 


3 DECT experiment at SSRF 


3.1 Beamline and CT system 


For X-ray imaging and biomedical application, the 
monochromatic X-ray CT imaging system in beamline 
BL13W1 of SSRF is shown in Fig.l. The two- 
(2D) CCD detector used 
experiment has 4000x2500 pixels, and its spatial 


dimensional in this 
resolution reached 13 um. Samples were placed on a 
rotation table approximately 30 m downstream source. 
Synchrotron radiation with a continuous energy 
spectrum was monochromatized with a double-crystal 
monochromator. BL13W1 is equipped with mono- 
chromator Si(111) and Si(311). In this experiment, 
Si(111) was used. The maximal photon flux density 
appears at the X-ray energy of 20 keV. A collimator in 
the system is used to prevent the scattering photon and 
improve the accuracy of reconstructed results. 

In DECT, it is intuitive to choose low/high 
energy as far from each other as possible. Because 


Table 1 Parameters of Beamline BL13W1 at SSRF 


Items 

Beam model 

Imaging method 

Energy range(KeV) 
Energy resolution(AE/E) 
Maximal beam size(mm) 


Photon flux density using Si(111) (phs/s/mm? @20keV@200mA) 


Max time resolution(ms/frame) 
3.2 Samples 


To evaluate the reconstructions and provide necessary 
correction, Water, alcohol, and polyethylene cylinder 
as phantoms are used to quantitatively evaluate the 
discrepancies between the measured and theoretical 
values. All the samples were scanned in plastic tubes 
of 18-mm diameter. 


too-low X-ray energy would increase image noise, the 
10 keV and 30 keV were chosen as the low/high 
energy. In order to ensure the same object position in 
scanning, each sample should take DECT twice under 
10 keV and 30 keV monochromatic X-rays alternately, 
and scan another sample. The low-energy X-ray 
absorption images were taken first within the m 
angular range under 10 keV. Then the sample was 
rotated backward to the beginning along the vertical 
axis. The X-ray energy was changed to 30 keV for 
next scanning. For each CT scan, the total number of 
views was 540 covering the m angle range. The data 
acquisition time and exposure, which depend on the 
photon flux density, takes about 2—4 min per CT 
scanning. It is necessary to slowly rotate the scanned 
object to prevent motion artifacts. The refraction effect 
can enhance the material boundaries and the image 
contrast, but also could cause inaccuracy in the 
electron density images. We shorten the object to 
detector distance at less than 10 cm to eliminate the 


boundary enhancement in the reconstruction images. 
b c = E g h ij | 
a | ; l | LF aar ee 
o | k 
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Fig.1 A schematic drawing of the experimental system for CT 
imaging in the beamline BL13W1 at SSRF. (a) White beam, (b) 
diaphragm, (c) filter 1, (d) y collimator, (e) white beam slit, (£) 
berrilium window, (g) filter 2, (hb) DCM (Double-crystal 
monochromator), (i) photon shutter, (j) y collimator, (k) rotation 
table, (1) 2-D X-ray detector. 


parameters 

Unfocused monochromatic beam 

PCI, Micro CT 

8-72.5 

<5x107 

45 (H)x5(V) 

1x10? 

1.0 

Water phantoms are used for calibration of 

electron density and correcting the systematic error. 
By water calibration, the reconstructed electron 
density value is multiplied by 1.03, which is the ratio 
of the measured to theoretical attenuation coefficients. 
The deviation mainly comes from the influence of 
noise and the nonlinearity in the detector response 
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function. Also, the theoretical process to derive from 


dual-energy linear attenuation coefficients may 
introduce error in the electron density measurement. In 
addition, the breast tissues in plastic tubes were used 
to investigate the feasibility of synchrotron radiation 


DECT in clinical application. 


3.3 Data processing 


Images without a sample are taken before and after 
each CT scanning to be Jp in reconstruction. X-ray 
energies of 10 keV and 30 keV are used in the DECT 
experiments. For each sample, two linear attenuation 
coefficient images with dual-energy X-ray are 
reconstructed by parallel-beam FBP method using the 
Shepp-Logan filter function after the dark-field and 
flat-field corrections. Electron density and effective 
atomic number of each pixel are calculated by the 
proposed reconstruction method. Solving Eq.(8), the 
coherent scattering term is calculated from standard 
element oxygen. The oxygen cross-sections were 
quoted from database”). 

In practical scanning, the data acquisition 
always suffers from noises and detector nonlinearity. 


Therefore, after retrieval of the attenuation coefficient 


maps from the raw projection data, total variation 
minimization method is used for noise reduction?" 
The computation from Eqs.(7) and (8) for the air part 
in the scanning FOV could result in a big unexpected 
error because the divisor is close to zero. Hence, we 
directly set the electron density in air region to a small 


value of 1x107. 


4 Results and Discussion 


In order to quantitatively evaluate the accuracy, we 
first used different cylinder phantoms including water, 
polyethylene and alcohol. Table 2 contrasts the 
The theoretical 


values of effective atomic numbers and electron 


measured and theoretical results. 


densities come from NIST standard reference database 
[20] For water and polyethylene phantom, the average 
value of measured electron density is fairly consistent 
with the theoretical ones (about 1% error). The 
electron density measured by synchrotron radiation 
DECT is sufficiently good, thus improving the 
precision of the treatment planning for radiotherapy. 
An effective atomic number image was obtained as a 
byproduct from DECT. 


Table 2 Experimental results of effective atomic numbers and electron densities of scanned phantoms 


Effective atomic number 


Electron density (g/cm?) 


Measured Theoretical Relative accuracy / % Measured Theoretical Relative accuracy / % 
Water 7.52 7.51 0.13 1.11 1.11 / 
Polyethylene 5.57 5:52 0.91 1.07 1.06 0.93 
Alcohol 6.53 6.47 0.93 0.90 0.89 1.12 


Figure 2 shows the electron density (a) and 
atomic number distribution images (b) of a breast 
tissue sample with pathological changes. The dark 
background is the air part in the tube. We can find a 
big difference between the two images. The image of 
electron density shows a striking contrast between 
normal and abnormal tissue, and are rather similar in 
the effective atomic number image. The information 
may give different criteria for the discrimination of 
abnormal tissues from normal tissues. Since the 
quantitative information about breast tissues is not 
currently available, we just calculated their average in 
different ROI. After statistics, the average electron 
density in the ROI area is 1.11 g/cm’; and 1.02 g/cm? 
in the ROI? area. 


Korz ©) 


T Plastic tube 
ROI 1 : 
a | Air 
` 


Fig.2 Electron density pe reconstruction results of breast 
tissue sample at the grayscale range of [0.8, 1.5] (g/cm°) (a). 
Effective atomic number Z reconstruction results of a breast 
tissue sample at the grayscale range of [1,15] (b). 
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5 Conclusions 


The precise measurement of electron densities is 
essential in treatment planning for hadron therapy. The 
most popular way is to convert the CT numbers to 
electron densities by establishing their relationship in 
spite of that it is not accurate enough. The SSRF has 
great potential for quantitative medical diagnosis and 
provide an excellent platform for X-ray imaging 
research. We conducted a DECT experiment to 
evaluate the accuracy and clinical application of 
with 
Different phantoms and 


electron density measurement synchrotron 


radiation. samples are 
investigated. The electron density reconstruction 
results of phantoms agree to the theoretical ones with 
high precision after correction, thus satisfying the 
heavy-ion and proton radiotherapy. Currently, the 
beamlines in SSRF is only available for research, and 
the SSRF should be designed for medical applications 
in the future. The DECT with synchrotron radiation 
has the potential to be an important tool in clinical 
performances. We will further study on this area, such 
as noise reduction methods, dose evaluation and 
clinical application research. 
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